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ABSTRACT 

Here we propose a mechanism for efficiently growing intermediate mass black holes 
(IMBH) in disks around supermassive black holes. Stellar mass objects can efficiently 
agglomerate when facilitated by the gas disk. Stars, compact objects and binaries can 
migrate, accrete and merge within disks around supermassive black holes. While dy- 
namical heating by cusp stars excites the velocity dispersion of nuclear cluster objects 
(NCOs) in the disk, gas in the disk damps NCO orbits. If gas damping dominates, 
NCOs remain in the disk with circularized orbits and large collision cross-sections. 
IMBH seeds can grow extremely rapidly by collisions with disk NCOs at low relative 
velocities, allowing for super-Eddington growth rates. Once an IMBH seed has cleared 
out its feeding zone of disk NCOs, growth of IMBH seeds can become dominated 
by gas accretion from the AGN disk. However, the IMBH can migrate in the disk 
and expand its feeding zone, permitting a super-Eddington accretion rate to continue. 
Growth of IMBH seeds via NCO collisions is enhanced by a pile-up of migrators. 

We highlight the remarkable parallel between the growth of IMBH in AGN disks 
with models of giant planet growth in protoplanetary disks. If an IMBH becomes 
massive enough it can open a gap in the AGN disk. IMBH migration in AGN disks 
may stall, allowing them to survive the end of the AGN phase and remain in galactic 
nuclei. Our proposed mechanisms should be more efficient at growing IMBH in AGN 
disks than the standard model of IMBH growth in stellar clusters. Dynamical heating 
of disk NCOs by cusp stars is transferred to the gas in a AGN disk helping to main- 
tain the outer disk against gravitational instability. Model predictions, observational 
constraints and implications are discussed in a companion paper (Paper II). 

Key words: galaxies: active - (stars:) binaries: close - planets-disc interactions- 
protoplanetary discs - emission: accretion 



1 INTRODUCTION 

Extensive evidence exists that supermassive black holes 
(> 10 6 Mq) are found in the c enters of most galaxies 
fe.g. iKormendv fc Richstondl 19951 ). Extensive evidence also 
exists for stellar mass black holes in our own Galaxy 
l|Remillard fc McClintockl 120061 ). Stellar mass black holes 
are expected to form as the end product of high-mass 
stars. Supermassive black holes, by contrast, have grown 
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to their current size over cosmic time, from much smaller 
■ 1 i 1 i— 1 



seeds (e.g. Begelman fc Reesl 1 19781: llslam. Taylor fc Silk.1 



|2004 iPortegies-Zwart et al.ll2004l ; iMicic et al.ll201ll , fc ref- 
erences therein). Intermediate mass black holes (IMBH; ~ 
10 2 — 10 4 Mq) may have been the original seeds for supermas- 
sive black holes or, they may have co ntributed to fast early 
growth of such seeds v ia mergers (e.g. iMadau fc Reesll200ll ; 
iMiller fc Colbertj|2004) . Though we expect IMBH should ex- 
ist, at least as an intermediate stage on the way to a super- 
massive black hole, observationally the evidence for their 
existence is scant and ambiguous, especially compared with 
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evidence for supermassive and stellar mass black holes. The 
low mass end of the supermassive black hole dis tribution in 
galac tic nuclei may extend down to ~ 10 s Mq ijjiang et al.l 
2011), but below this mass the evidence becomes ambigu- 
ous. The ultra- luminous X-ray sources (UL Xs) observed out- 
side galactic nuclei (e.g. IWinter et a l. 2009) may b e powered 
by accretion onto IMBH l|Miller fc Colbert! |2004 ). However 
ULXs could also be a explained b y beamed ra diation from 
accreting stellar-mass black holes (King 2009) and power- 
law dominated ULXs might be due to background AGN. 
IMBH have so far been hard to find and constrain in the 
local Universe, either in our own Galaxy or at low z. 

Active galactic nuclei (AGN) are believed to be powered 
by accretion onto a supermassive black hole. The accretion 
disk should contain a population of stars and compact ob- 
jects (collectively nuclear cluster objects, NCOs) that can 
migrate within and accrete from the disk (e.g. Ostrikerl 
19831; ISver. Clarke fe Reesl Il99ll; lArtvmowicz et all Il993al; 
Goodman fc Tan! |2004|; iLevinl 1200?); iNavakshin fc Sunvaeyl 
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20071 ; iMcKernan et alj|2011al lri). In lMcKernan et all (|2011a 
we speculated that IMBH seeds may form efficiently in AGN 
disks due to NCO collisions and mergers, which is quite dif- 
ferent from the standard model o f stellar mass black holes 
merging in stellar clu sters (e.g. iMiller fc Hamilton] |2002| ; 
IMiller fc Colbertll2004h . Here we argue that IMBH produc- 
tion is in fact far more likely and more efficient in AGN 
disks, with implications for AGN observations, duty cycle 
and supermassive black hole accretion rates. 

In this paper (and its companion, Paper II, McKernan 
et al. 2012) we discuss semi-analytically the production of 
intermediate mass black holes in the environment of AGN 
disks. Discussion of observational predictions of this model 
of IMBH growth as well as consequences for AGN disks, duty 
cycles and the demographics of activity in galactic nuclei at 
low and high redshift will be left to Paper II. In section fj2j 
we discuss why we think IMBH can be built in AGN disks. 
In section fJ3]we explore mechanisms that will be important 
in actually growing IMBH in AGN disks, including the com- 
peting forces of eccentricity damping and excitation in the 
disk. The importance of IMBH migration is outlined in sec- 
tion fj4] Section ij5] outlines a simple model of IMBH growth 
in AGN disks and we highlight the remarkable parallel be- 
tween the growth of IMBH in AGN disks and the growth of 
giant planets in protoplanetary disks. Finally in section SjfJ 
we outline our conclusions and future work. 



2 WHY IMBH CAN BE BUILT IN AGN DISKS 

The largest, supermassive, black holes in the Universe 
(Mbh ~ 10 s — 10 9 M Q ) l ive in galactic centers (e.g. 
iKormendv fc Richstond [19951. We expect a dense nuclear 
cluster of objects to surround the supermassive black 
hole as a result of stellar evolution, dyna mical fr i ction , 
secular evolution and minor mergers (e.g. iMorris! 1 19931 ; 
iMiralda-Escude fc Gouldl |2000| ; iMerrittl [201ol ). In our own 
Galaxy, the distributed mass within ~ lpc of Sgr A* is 
~ 10 — 30% of the mass of the supermassive black hole 
|Schodel et all 120091 ) . If a large quantity of gas somehow 
arrives in t h e innermost pc of a galactic nucleus (e.g . 
Krolikl 1 19991; iKawawatu et al.l [20031; IVittorini et al.l 120051; 
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likely lose angular momentum and accrete onto the central 
supermassive black hole. But in doing so, gas must also in- 
teract with the NCO population. Depending on the aspect 
ratio of the disk that forms, a few percent of NCO orbits are 
likely to coincide with the accretion flow. A small percentage 
of NCO orbits coincident with the geometric cross-section of 
a thin disk would lead to an initial population of ~ 1O 3 M0 
of NCOs in a pc-scale accretion flow around a SgrA* sized 
black hole. 

NCOs can exchange angular momentum with gas in the 
disk, and each other, so t hey can scatter each oth er and mi- 
grate within the disk fsee lMcKernan et al.ll20l"lah . The pro- 
cesse s involved are analagous to protopla netary disk theory 
(e.g. iPollack et al]|l996l ; TAr"mitagell2010l ). Indeed, physical 
conditions in the outskirts of AG N disks are relatively clos e 
to those in protoplanetary disks (jMcKernan et all l2011aT ). 
The migration of NCOs in the disk will enhance the proba- 
bility of collisions, mergers and ejections. Under these con- 
ditions, IMBH seeds can grow. IMBH seeds will be objects 
> IOMq that will not lose very much mass (e.g. stellar mass 
black holes, hard massive binaries). IMBH 'seedlings' we de- 
fine as objects ^ IOMq that have grown via mergers (e.g. 
the merged end-product of a hard binary). IMBH seeds and 
seedlings located at semi-major axis a in an AGN disk will 
maintain a 'feeding zone' within which they may collide with 
nearly co-orbital disk NCOs. By analogy with proto-planet 
growth we define the feeding zone to be a ± 4Rh where a is 
the IMBH semi-major axis and R H = a(g/3) 1/3 is the IMBH 
Hill radius, with q the mass ratio of IMBH:supermassive 
black hole. Once the object gets to ^ lOOM© we will call 
the result an IMBH. 

The disk NCO population is subjected to dynamical 
heating from cusp stars and dynamical cooling from gas 
damping. If gas damping dominates, IMBH seeds and disk 
NCOs will have their orbits rapidly damped. As a result 
their collision cross-sections will rapidly increase (since the 
relative velocity of encounters will be small, particularly in 
the outer disk). IMBH seedlings will initially accrete disk 
NCOs within their feeding zone in a 'core accretion' mode 
of growth. Once nearby NCOs have been scattered or ac- 
creted, gas accretion dominates IMBH growth. However, the 
migration of IMBH seedlings within the disk allows growth 
to continue via collisions as well as via gas accretion. Thus, 
we expect IMBH to grow within AGN disks, analagous to 
the growth of giant planets within protoplanetary disks and 
we expect the IMBH growth rate will be much larger than 
in stellar clusters. 

Here we concentrate on growing IMBH within the AGN 
disk itself. Of course it is possible that IMBH already exist 
in the galactic nucleus when the AGN disk first forms. A 
top-heavy initial mass function of cusp stars can lead to 
IMBH seedling formation before low angular momentum gas 
arrives in the nucleus. IMBH can also arrive from outside the 
galactic nucleus to interact with the AGN disk since mass 
segregation and dynamical friction can deliver IMBH to the 
central parse c of galactic nuclei in a few Gyrs from nearby 
clusters (e.g. IMcKernan et al.ll2011bl . see also Paper II). 

The physics involved in IMBH formation in AGN 
disks spans multiple regimes and physical processes and 
would usefully benefit from detailed numerical simula- 
tions. Such simulations require realistic treatments of 
(amongst others): N-body collisions, mergers, accretion, 
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tidal forces, gravitational radiation, Special & General Rel- 
ativity, radiative transfer, the magneto-rotational instabil- 
ity and the gravitational instability. However at present 
there are no simulations that can adequately address the 
relevant physics in a self-consistent manner. We take a 
semi- analytic approach following (e.g. iMiller fc Hamilton! 
20021) on the build-up of IMBH in star clusters, (e.g. 
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Alexander. Begelman fc Armitagel |2007| ) on stellar dynam- 
ical heating and cooling, as well as f ormalism on planet 
growth from pr otoplanetary theory fe.g. IPollack et aTJl 19961 : 
lArmitaga l2010. & references therein). 



3 HOW TO BUILD IMBH IN AGN DISKS 

In this section, we shall outline the key phenomena involved 
in growing IMBH seeds in AGN disks. In order to grow 
into IMBH, seeds must collide with and accrete mass, ei- 
ther NCOs or gas. In section H3.ll we discuss NCO collision 
cross sections in AGN disks and the importance of eccen- 
tricity damping and excitation. In section H3.2I we outline a 
model of dynamical heating and cooling of NCO orbits in 
AGN disks and we discuss the implications for IMBH seed 
growth. Section H3.3l brieflv outlines issues involved in merg- 
ing binaries, a potentially important channel for producing 
IMBH seedlings. 



3.1 Collision cross-sections in the disk 

In the absence of disk gas, NCOs change their orbits only due 
to weak gravitational interactions, occuring on the (long) 
relaxation timescale (see below). The interaction with a 
gaseous disk gives rise to new effects, namely torques. NCOs 
embedded in the disk, or crossing it, will have their eccen- 
tricities and inclinations damped relative to the disk. Such 
processes should enhance the stellar density in the disk re- 
gion and lower the velocity dispersion of NCOs embedded 
in the disk, in the absence of other important effects. This, 
in turn, gives rise to a higher rate of encounter and colli- 
sions between NCOs in the disk. The collisional cross-section 
(o"coii) of compact NCOs of mass M depends on the relative 
velocity at infinity (i>oo) as 



croon ~ nr p (2GM/vl ) 



(1) 



in the gravitational focussing regime, where r p is the sep- 
aration at periastron. In AGN disks the relative velocities 
involved in close encounters can be very small compared to 
the velocity dispersion in star clusters (typically ~ 50km/s). 
For example, NCOs on circularized orbits separated in the 
disk by AR ~ 0.01-R at R = 10 5 r 9 have relative veloci- 
ties due to Keplerian shear at periastron of only ~ 5kms~ , 
where r g — GMsmbh/c 2 is the gravitational radius of the 
supermassive black hole of mass Msmbh- The disk NCO 

velocity dispersion varies with radius as a ~ \/e 2 + ~?Vk, 
where e, i, Vk are the mean NCO orbital eccentricity, mean 
NCO orbital inclination and Keplerian velocity respectively. 
Fig. [1] shows the NCO velocity dispersion (a) as a func- 
tion of radius in a Keplerian AGN disk for a range of mean 
eccentricities and inclinations. Also shown in Fig. \T\ is the 
typical velocity dispersion in star clusters (~ 50km s , red 
horizontal dashed line). So Voo for a typical interaction in 




10 4 10 s 
Radial Distance (Rj 



Figure 1. The velocity dispersion (cr ss \J e 1 + i 2 Vfc ep ) of NCOs 
in a Keplerian AGN disk. Shown are a as a function of disk 
radius, for eccentricity and inclination values of (e,i)=0.01,0.05 
(solid lines) and e = 0.01, 0.05, 0.1, with i = 0(dashed lines). Also 
shown (red dashed horizontal line) is the typical velocity disper- 
sion in star clusters (~ 50 km s" 1 ). Note that most disk NCOs 
should live in the outer disk (> 10 4 r a ) for an NCO population 
that grows as r 2 . 



a stellar cluster is ~ 50km s^ 1 . Fig. [T] shows that the ve- 
locity dispersion of NCOs at large disk radii is less than in 
star clusters for small to moderate NCO orbital eccentrici- 
ties and inclinations (e, i ~ 0.01 — 0.05). Since the numbers 
of NCOs should increase with radius, most NCOs should live 
in the outer disk, where the NCO velocity dispersion should 
be smallest. 

The collisional cross-section of a seed IMBH (mass M) 
with compact objects (mass m) such as neutron stars, white 
dwarfs a nd stellar mass black hole s depends on relative ve- 
locity as (|Quinlan fc Shapi"rol ll989) 



O"coll = 2-7T 
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Numerically, this 
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be written 
i ln >;-„ ' ,; liii where Uoo = 10 B W6cms 
M5Q,mio are in units of 50Mq , IOMq respectively 
i|Miller fc Hamilton! [2002] '). For M 50 ,m 10 = 1, located ~ 
10 5 r 9 from a supermassive black hole, small eccentricities 
e ~ 0.01 in Fig. [T] lead to collision cross-sections up to an 
order of magnitude larger than in clusters (for Mso,mio = 
1;«6 = 5 above). However, for large NCO orbital eccen- 
tricities (e ^ 0.1), IMBH collisions in AGN disks will have 
smaller cross-sections than in star clusters, over most of the 
disk. Therefore if the initial mean eccentricity (e) of the disk 
NCO distribution is large, mechanisms for damping orbital 
eccentricities will be very important in determining whether 
IMBH growth via collisions is efficient in AGN disks. 
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3.2 NCO orbital damping & excitation 

We begin with a fully analytic approach, demonstrating the 
relative importance of competing terms and effects. The ve- 
locity dispersion (a) of NCOs in a disk is excited by dynam- 
ical heating and is damped by dynamical cooling. Thus 

^ = AQ+ - AQ- (3) 

where AQ+ is the dynamical heating term and AQ- is the 
dynamical cooling term. Dynamical heating comes from two 
sources: the relaxation of disk NCOs through mutual inter- 
actions and the dynamical excitation of disk NCOs by cusp 
NCOs, so 



AQ + = <5Q r clax + ' 



(4) 



Considering first the relaxation term, we assume that there 
are Ni stars of mass Mi and velocity dispersion <xi in an an- 
nulus of wi dth AR centered on R. The relaxation ti mescale 
is given by l|Alexander. Begelma n fc Armit agc 2007) 



^relax — 



5Qr 



2TrCiRARa i 
G 2 Nm 2 lnAn 

a _ Lh J_ 

^rclax Cl <T 3 



where 



G 2 iViMi lnAi 



(5) 
(6) 

(-) 



where fl is the Keplerian frequency, lnAi (~ 9) is the 
Coulomb logarithm and C\ ~ 2.2. The solid curve in Fig. [2] 



shows the evolution of < e 
for a population of TVi = 10' 
in an annulus of width AR = 



1/2 



due to relaxation alone 



stars of mass M\ — 
O.lpc centered on R 



O.6M 
= O.lpc 



(equivalently 1 - 3 x 10 r g of an AGN disk around a 10 8 M Q 
supermassive black hole). For these stars, Vk = 2100km/s 
and so t or b = 9 x 10 9 s. Since, for moderate eccentricities, 



< e 2 > 1 / 2 v k 



(8) 



the solid curve in Fig. [2] follows a t 1 ^ form and rises (limited 
by the increase in a to approximately the Keplerian veloc- 
ity). This solid curve applies to an isolated annulus of stars 
in the absence of competing effects. 

Additional heating is supplied by the cusp population. 
The cusp stars will transfer kinetic energy to the 'colder' 
disk popul ation a n d exc ite the o\ distribution of the disk 
NCOs Cse elPeretsI <|200S|) and Perets et al., 2012, i n prep.). 
Following lAlexander. Begelman k, Armitagel l|2007l ) 5Q cxcitc 
has the form 



scxcitc — 



c 2 



(9) 



where the cusp population Ni 2> Ni has a similar mass 
function Mi = Mi to the NCOs in the disk and 



Dz 



G NjMjMilnAi 

RiARt or b 



(10) 



with <7ii = (cti + (ii)/2 and Eij = SMi^o"^ is the kinetic 
energy. Note that £ eX cite is analagous to the i rc iax term in 
equation [5] except is now the cusp population (Ai).So 




Figure 2. The rms eccentricity as a function of time for an annu- 
lus of 10 3 stars of identical mass 0.6Mq located in an AGN disk 
between 0.05-0.15pc (or 1 — 3x 10 4 r 9 ) around a 10 s Mq supermas- 
sive black hole. < e 2 > x / 2 scales as TV 1 / 4 , M 1 / 2 . The solid curve 
shows the relaxation of the isolated distribution of stars assum- 
ing initial circularized orbits. The dashed curve shows relaxation 
plus dynamical heating from a population of Ni = lOOTVi cusp 
stars with Mi = Mi and "&i ~ 0.5iifc. The dotted line shows the 
net domination of exponential damping by disk gas over orbital 
excitation by cusp stars. 



toxcito ~ (Ni/Nijtrei^ < t rc iax- We assume that on aver- 
age cfi ~ \/e 2 + iVk ~ 0.5vk so Ei > Ei (i.e. the cusp 
stars have greater kinetic energy than the disk stars). The 
dashed curve in Fig. [2] shows the addition of this excitation 
term to the evolu tion of < e 2 > 1 ' 2 , assuming 
and Ci/Ca = 3.5 j Alexander. Begelman fc Armitagell2007l ). 
with InAii ~ lnAi. Clearly, the curve retains a t 1 / 4 depen- 
dence, but at larger values of < e 2 > 1//2 . Thus, dynamical 
heating by stars in the cusp domina tes relaxation by the 
stars in the disk (see also |Peretsll2008l ). 

The competing dynamical cooling term AQ- is dom- 
inated by gas damping of the disk NCO orbits. Gas drag 
in AGN disks will tend to reduce small NCO orbital ec- 
centricities and inclinations to much smaller values. Gas 
at co-orbital Lindb lad resonances will damp (e,i) for NCOs 
with q $ 10~ 3 (e.g.lArtvmowicz et al.lll993bl;IWard fc Hahnl 
11994 ICresswell et alj|2007l : iBitsch fc Klevll201Ch . Since this 



mechanism depends on the co-rotating gas mass, both stel- 
lar and compact NCOs and IMBH seeds will have their 
orbits damped, particularly in the outer disk where most 
of the disk mass is located. For small eccentricities (e < 
2(H/r)), orbital eccentricity decays exponentially over time 



(H/r) 



where h = H/r is the disk aspect ratio 



and T m j B is the mig ration timescale l|Ward Sz H ahn 1994; 
IBitsch fc Klevll201Ch . Thus 



de 

~dt 



(ii) 
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which will give us a term linear in a as the damping term in 
AQ-. By analogy with the relaxation term (cr/irelax) above, 
we c hoose k = l/tdam p- The damping timescale is given by 
(e.g. iHorn et alj|20l3 ) 



M 2 

mEa 2 S7 



(12) 



where h gas = H/R — c s /vk is the gas disk aspect ratio (the 
disk of stellar NCOs has an aspect ratio ftstars = cr/vk but 
the disk NCOs are not the source of damping). For larger 
eccentricities (e > 2( H/r)) the eccentricity damping goes as 
i|Bitsch fc Klevll2O10D 



K 



dt e 

where we choose the same normalization n — 1/tdamp 
above. So, 



(13) 



AQ- = 



ft 



(14) 



where /?' = 1 if e < 0.1, zero otherwise and /?" = 1 if e > 0.1, 
zero otherwise. Thus, our expression for the combined relax- 
ation, excitation and gas damping of the velocity dispersion 
for an annulus of NCOs is given by 



da 



Di J_ D2_oj_ f,Ei\_ 
d af + C 2 a\ t \ Ej K 



0a- 



P" 



(15) 



For < e 2 > 1/2 < 2h (|Bitsch fc Klevll2010l ), eqn. [15] has the 
general form 



da 
~dt 



^rcla 



A 



tda 



= 77 — KO". 



(16) 



where we have combined the dynamical heating terms into 
a single general A /a 3 term. Since t CX citc <S trciax, at steady 
state, da/dt = 0, and so tdamp ~ icxcitc and e takes the 
general form 



4G 2 N c mM£ u /ig as lnA 
27vC 2 T,a 2 RARvi 



(17) 



where jV c is the number of stars in the cusp. So, for a popula- 
tion of 10 3 x 0.6M Q stars located at 1 - 3 x 10 4 r 9 in an AGN 
disk, with /i gaa ~ 10 -2 and a cusp population of N c — 10 
stars, equilibrium eccentricity is e ~ 0.01. In Fig. [5]we plot 
(dotted line) the evolution of e over time assuming i ~ 0. 
From this we see that disk NCOs should rapidly settle down 
to near circular orbits (< e 2 > 1//2 ~ 0.01) within ~ O.lMyr. 
Therefore collision cross-sections (cr C oii) of IMBH seeds in 
AGN disks should rapidly become much larger than typical 
collision cross sections in star clusters and it is gas damping 
that makes the difference. 

We expect gas damping to become even more dominant 
if the NCO disk population declines (N-) due to mergers, 
accretion and scatterings. From equipartition of energy, we 
expect N- will be dominated by low mass stars. At moderate 
inclinations, these NCOs can be captured fairly quickly by 
the disk again (such that N+ in creases), if the gas damp s 
the orbital inclination efficiently IjArtvmowicz et al.lll993al ). 
As N+ increases, the system is driven towards a dynamical 
equilibrium when N- ~ N+. 

One important point to note from the above discussion 
is that the dynamical heating of the NCOs by cusp stars 



(AQ+) gets transferred to the AGN disk gas. The stabil- 
ity of the outskirts of the AGN disk is a well -known and 
unsolved problem (e.g. lSirko fc Goodman! [20031 ); dynamical 
heating of disk NCOs by cusp stars is a new, additional 
source of disk heating which will contribute to maintain- 
ing the outer disk against gravitational instability. A self- 
consistent calculation of the disk heating requires a disk 
model (e.g. lSirko fc Goodmanll2003h and is beyond the scope 
of this paper, but see McKernan et al. 2012 (in prep.). Never- 
theless, we can see that a large density of NCOs in a galactic 
nucleus will strongly excite the orbits of disk NCOs (<5Q e xcite 
is large). Gas damping (AQ-) will naturally transfer much 
of this dynamical energy to the disk gas. Thus, disk luminos- 
ity must increase and the disk itself will puff up. The scale 
height increase will be a function of the density of NCOs in 
the nucleus. Therefore, one prediction of our model is that 
among nuclei with similar supermassive black hole masses, 
those with denser stellar cusps, should generate more lumi- 
nous AGN disks (see Paper II). 

So far we have discussed low mass stars. However, we are 
interested in higher mass IMBH seeds. For simplicity let us 
assume a steep NCO mass function (dN/dM oc M~ 3 ) with 
two mass bins. The low mass population NCOs are 0.6Mq 
stars (Ni in number); thus the high mass NCO population is 
10 _3 7V; xlOM© stellar mass black holes. For a total initial 
disk NCO mass of 10 3 (10 4 )M Q , the distribution is 1.65 x 
10 3 (10 4 ) low mass stars and 1(10) stellar mass black holes. 
lAlexander. Begelman fc Armitagg (|2007l ) show that a low 
mass population of stars will diffuse out of the disk more 
than the high mass population of stars and in fact damp the 
orbits of the high mass stars, as expected from equipartition. 
Thus, for small initial values of < e 2 > x / 2 among disk NCOs, 
we expect potential IMBH seeds in AGN disks to evolve 
to even smaller eccentricities than the equilibrium value of 
e ~ 0.01 calculated above. Recall that small eccentricities 
imply large (j co ii, allowing IMBH seedlings to grow rapidly 
via collisions. 



3.3 Binary mergers in the disk 

Depending on the recent star formation history of a given 
galactic nucleus, massive binaries are likely to be rare in 
AGN disks. However, if there is even one in the initial AGN 
disk, it will have the largest collisional cross-section of any 
disk NCO and should undergo the largest number of interac - 
tions (e.g. iPortegies-Zwart et al.lfl999l : iFregeau et al.|[2004l ). 
A massive binary, if present, is therefore the most likely 
IMBH seed and should arise frequently enough to be of as- 
trophysical interest (for similar ideas concerning planetes- 
imal growth through b inary-single in teractions in a proto- 
planetary disk see e.g. iPeretj l|201ll )). In this section, we 
briefly consider some of the issues involved in binary merg- 
ers in an AGN disk and we contrast the merger efficiency 
with that found in star clusters. For eas e of comparison we 
consi der the 50M© + 1OM binary from dMiller fc Hamiltonl 
120021) . 

An unequal mass binary (M > m, separation at>i n ) in 
the disk is considered hard if its binding energy (GMm/abin) 



K) 



of 



is greater than the kinetic energy (ma = m(e + % 
a typical interacting NC O. The collis ional cross-section of 
such a binary is given by (|Peretsl l201ll) 



6 B. McKernan, K.E.S. Ford, W.Lyra & H.B. Perets 



ffcoll 



2 

Tia b 



4-)' 

\Voo / 



r N /10 6 



Obin/2.14 x 10 8 cm 



(18) 



where v c = (G//x(Mm/a b i n )) is the critical velocity sepa- 
rating hard and soft binaries, with /j, = (Mbm x Mn)/ (Mbi n + 
Mn) is the reduced total mass of the binary (Mbm) and in- 
teracting NCO (Mn). The time to merge for a binary of re- 
duced mass [i = mM/(M + m) where M > m, semi-major 
axis Obin and eccentricity e bin is 

r mcrgc « 3 x W 8 Ml ) (tJ,M 2 )- 1 (a bi n/Re) 4 (l - e 2 bin ) 7/2 yr (19) 

and the typical semi-major axis separation for a merger time 
of TeMyr (assuming ebin ~ 0) is 

a bin ~ 3 x lQ U T 6 1/4 M^ 2 m^ 4 cm (20) 

where Mso, mio are the masses in unit s of 50Mq and lOAf© 
respectively (jMiller fc Hamilton! 12002 ). However, the above 
disc ussion neglects the gas disk. 

iBaruteau et al.l (|201ll ) carried out hydrodynamic simu- 
lations of a binary in an AGN disk and found that binaries 
harden rapid ly due to interac t ion w ith their own migratory 
spiral wakes. IBaruteau et al.l (|201lF ) also found that ebin is 
damped rapidly with inward migration. The rate of binary 
hardening (dbin) scales with the disk surface density such 
that abin/obln "C Tbin, the binary migration timescale. Mas- 
sive binaries with initial separation abin ~ O.SRh, end up 
at half this separation within 10 orbits of the supermassive 
black hole, where Rh is the Hill radius = (g/3) 1//3 a, with 
a the semi-major axis of the binary center of mass and q 
the mass ratio of the reduced mass binary to the supermas- 
sive black hole. For a constant rate of hardening (& b i n ), a 
M50, mio migrating binary separated by ~ 3 x 10 11 (10 13 )cm 
(or 3i? Q (2 AU)) at 10 5 r 9 will merge in < 0.1(2)Myrs, which 
is a very small fraction of the AGN disk lifetime. So, because 
of the presence of a gas disk, it is easier to harden binaries 
in AGN disks than in star clusters. 

Binaries will also encounter field NCOs in the disk as 
they migrate. This can either result in hardening to merger 
or disruption. The probability that a binary is disrupted per 
unit time is 1 /tdis where 



^dis — 



9\E\ : 



16^vG 2 m 4 a 



1 + 



4mcr 



1 + exp 



( Z\E\ 
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(21) 



where the field NCO has mean mass m and velocity dis- 
persion a, v is the number density of field NCOs, E — 
—GMm/aun is the binding energy of the binary, and 
the average energy chang e per interaction is ~ — 0.2m<r 2 
ijBinnev fc Tremainelll987l l. In stellar clusters ~ 10 2 field in- 
teractions are required to harden massive binaries to merger 
jMiller fc Hamilton! [2002) , where the relative velocities at 
close encounters are approximately the velocity dispersion 
(~ 50kms _1 ) in star clusters. In AGN disks, the number of 
interactions required to harden a binary to merger depends 
on the mean eccentricity (e) of the field NCO orbits. The 
probability of binary disruption (1/idis) increases as e in- 
creases. However, for already hard binaries (large \E\), fewer 
interactions (of energy — 0.2mcr 2 ) are required for merger in 
gas disks. Thus, hard binari es will continue to ha rden due 
to inward Type I migration (IBaruteau et al.ll201ll ) and due 
to interactions with field NCOs in the disk. 

For a 6OAf unequal mass binary (Mso,mio) at 10 5 r 9 
interacting with a population of NCOs having e = 0.01, the 



energy per binary interaction (— 0.2mcr 2 ) is ~ 4% the typical 
interaction energy in stellar clusters (where a ~ 50km s _1 ). 
So for highly damped NCO orbits, binary interactions are 
likely to be 'soft'. The number of interactions depends on 
the disk NCO surface density. If this binary is already hard, 
it could merge within a few orbits, i.e. orders of magnitude 
faster than binary merger timescales in stellar clusters. 



4 NCO MIGRATION & COLLISIONS 

The gas in the AGN disk exerts a net torque on NCOs. 
This means that individual NCOs (and IMBH) will migrate 
within the gas disk, enhancing the probability of collision 
and merger. Migration is mostly inward in disks, although 
sometimes outward. This means that IMBH can migrate into 
new regions of the disk, in search of new NCO 'victims'. The 
IMBH feeding zone (approximately a ± ARh) moves with 
the migrating IMBH. This is analagous to a giant planet 
core continuing to collide with planetesimals as it migrates 
through a protoplanetary disk lAlibert et all l|2004h . Migra- 
tion can stall in disks, leading to a pile- up (overdensity of 
disk NCOs). In this case rapid merger leading to rapid IMBH 
growth ma y occur, by analo gy with pile-up in protoplane- 
tary disks (|Horn et al.ll2012l ). A detailed calculation of this 
scenario will be carried out in future work. NCOs and IMBH 
may also migrate onto the central supermassive black hole, 
just as protoplanets may migrate onto a central star. So the 
issue of IMBH survival in AGN disks parallels the survival 
of giant planets in protoplanetary disks. 



4.1 Type I IMBH migration 

NCOs with mass ratios q ^ 10 -4 of the mass of the cen- 
tral supermassive black hole will undergo migration in the 
disk analago us to Type I protoplanet ary migration, on a 



disk analagous to type 1 protoplaneta 
timescale of |Paardekooper et al.ll201oh 



n 



1 M 

N qT,r 2 



\ r J oj 



(22) 



where M is the central mass, g is the ratio of the satel- 
lite (NCO) mass to the central (supermassive black hole) 
mass, E is the disk surface density, H/r is the disk aspect 
ratio and u is the satellite angular frequency. The numer- 
ical factor N depends on the ratio of radiative to dynam- 
ical timescales and is a function of the power-law indices 
of S, T and entropy (|Lvra et al.l |2010| ; iPaardekooper et al.l 
l2010f) . Note that the Type I migration timescale decreases 
with increasing migrator mass at a given radius so more 
massive NCOs will migrate more quickly at a given disk ra- 
dius. Binary NCOs face exactly the same torques and will 
migrate on the same timescale, but q and r in eqn. [22] are 
replaced with, the ratio of the reduced mass of the binary 
to the supermassive black hole and a, the location of the 
bin ary center of mass i n the A GN disk respectively. 

ISirko fc Goodman! (|2003h model an AGN disk includ- 
ing all the parameters we require for calculation of migrator 
timescales as a function of radius aroun d a 10 8 M(T) super- 
massiv e black hole. Although in principle lSirko fc Goodman! 
(2003) model a disk out to 10 7 r 9 , they regard their disk as 
effectively truncating at ~ 10 5 r 9 . This disk also requires a 
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constant mass accretion rate (M) and a constant disk viscos- 
ity (a) at all radii over the disk lifetime, which are obvious 
simplifications. Neverthe less using the simple AGN disk of 
ISirko fc Goodman! (|2003) as our disk model, we can estimate 
migration timescales semi-analytically. Figure [3] shows the 
Type I migration timescales of a fiducial IMq NCO (upper 
curve) and a 60Mq IMBH seed (lower curve) as a function 
of disk location. The curves in Fig.[3]are generated by choos- 
ing N ~ 3 in eqn !22l and assuming that S and H/r have th e 
form of the curves in Fig. 2 of Isirko fc Goodman! (|2003h . 
Also marked in Fig. [3] is an approximate AGN lifetime of 
50Myrs (red dashed line). Evidently, substantial changes of 
NCO orbital radius can occur even for low mass NCOs in the 
inner disk (< 10 4 r 9 ) over the AGN lifetime (few xlOMyr). 
Larger mass migrators (stellar mass black holes, binaries, 
large mass stars or seed IMBH) are likely to have migration 
timescales roughly comparable with the AGN disk lifetime. 
Therefore, as IMBH seedlings grow in mass they should mi- 
grate in the disk and encounter low mass NCOs at low rel- 
ative velocities. If migration stalls for a given NCO, either 
inwards at small disk radii, or outwards at large disk radii, 
interactions are possible at incredibly low relative velocities 
(voo). Note that spiral density waves from migrating NCOs 
should not strongly per turb the NCO migrations or their 
orbits (|Horn et alj|2012l) . Here we assume that seed IMBH 
migrate independently. Of course, resonant capture can oc- 
cur, both between IMBHs and NCOs and between multi- 
ple IMBHs, analagous to resonances between Jupiter and 
its moons, or between Neptune and Pluto. Given the low 
mass ratios and high migration speeds, an assumption of 
independent migration seems reasonable, but future simula- 
tions involving multiple migrators are required to test this 
assumption. 

Because we expect NCOs in the AGN disk to migrate 
differentially, we expect the migrators to encounter each 
other. As large mass NCOs migrate inwards across the or- 
bits of less-massive NCOs, if the gas has damped e (see 
SJ3TT] above) the relative velocities («oo) will be low and the 
collision cross-section (with gravitational focussing) will be 
large. However, different large mass NCOs will have differ- 
ent outcomes from multiple interactions. Large stars could 
have a shorter life expectancy if there are many NCO inter- 
actions increasing the odds of merger and supernova. Stellar 
mass black holes and IMBH seedlings will tend to shred low 
mass main sequence or giant stellar NCOs as they migrate 
inwards, although they can swallow compact NCOs whole 
as tidal forces shred compact objects only after the compact 
object crosses the innermost stable circular orbit (ISCO). 
Inward migrating binaries will tend to harden and scatter 
lower- mass NCOs until merger (see below). 

Although a majority of Type I migration is directed 
inwards, Type I migration can also occur outwards in proto- 
planetary disks. For eccentricitie s e < 0.02, migration may 
be outwards rather than inwards jBitsch fcK lcv 2010). Mi- 
gration can also stall on both inward and outward migra- 
tions, depending on the temperature and density of the ad i- 
abatic disk (e.g. IVeras fc Armitagell2004 iLvra et alj|2010h . 
So, as we consider a population of NCOs migrating in an 
AGN disk and interacting with each other, we do so with 
the caveat that a fraction of the NCOs may be migrat- 
ing outwards, or stalled. From equation 1221 the migration 
timescale gets longer at small disk radii (r), where £ de- 
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Figure 3. Estimates of Type I migration timescales for IMq 
NCOs (upper curve) and 60Mq IMBH seedlings (lower curve) as a 
functi on of radius in the AGN disk modelled bv lSirko fc Goodman! 
(2003). Also shown (red dashed horizontal line) is a fiducial AGN 
lifetime of 50Myrs. Most disk NCOs should live in the outer disk 
(> 10 4 r 9 ) for an NCO population that grows as r 2 . Note that 
migration is fastest around 10 3 r 9 in the disk, where E is largest 
and the disk is thinnest (H/r is smallest). There is a possibility 
of NCO pile-up due to migration at smaller disk radii, as the 
co-rotating mass of gas drops dramatically in the disk interior. 



creases and (H/r) increases due to disk heating. Migration 
may even stall or cease at small disk radii, particularly for 
large migrator masses, when the co-ro tating disk mass be- 
comes less than the migrator mass (e.g. ISver Sz Clarkdll995l : 
lArmitagell2007h . However, conditions in the hot inner disk 
may be dramatically different to the outer disk. If there is 
an abrupt transition to an optically thin accretion region, 
or a disk trucation or cavity, migration will stall. In this 
case, NCO pile-up can occur, potentially leading to mergers 
and ejections. Recent N-body simulations of protoplanetary 
disks suggest that migrator pile-up in regions of the disk 
where inward and outward torq ues balance results in very 
rapid merging (|Horn et alj|2013 ). If migrator pile-up occurs 
in AGN disks, it could favour the rapid building of IMBH 
seedlings. A stalled IMBH seed can merge with and scat- 
ter piled-up NCO migrators. Using a simple equipartition of 
energy, an IMBH seedling of mass M stalled at 10 2 r 9 in an 
AGN disk with e, i ~ 0.01 could scatter in-migrating NCOs 
(m) to a ~ (M/m)400 km s" 1 . At such hypervelocities, 
small mass NCOs can be ejected into a galactic halo (see 
Paper II). 

If disk NCOs migrate inwards, their rate of migration 
decreases in the inner disk (< 10 3 r 9 ) as the disk surface den- 
sity and co-rotating disk mass drops (see also Fig. [3] above). 
Although conditions in the innermost AGN disk are dra- 
matically different from those in the outskirts, it is useful 
to consider the possibility of NCO pile-ups. Evidently, if the 
AGN inner disk truncates at some radius and then becomes 
a geometrically thick, optically thin advection flow, NCO 
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migrators will stall near the disk truncation radius. As NCOs 
build up over time, the chances of interaction increase and 
hypervelocity scatterings become likely. The conditions may 
be similar to the migration trap for N protoplan ets in proto- 
plane tary disks, where rapid merger is possible l|Horn et al.l 
12012ft . In this case, IMBH could form in the inner disk with 
distinctive observational signatures (see Paper II for details) . 
Of course, migration traps can also occur as a result of out- 
migration. 

A majority of NCOs in galactic nuclei will have or- 
bits that do not coincide with the plane of the AGN disk 
(i ^ 0.05) and will instead punch through the disk pe- 
riodically. These NCOs can interact with each other, the 
disk and the migrating NCOs in the disk. NCOs on or- 
bits with small radii will eventually de cay into the plane 
of the disk over the AGN disk lifetime (|Artvmowicz et akl 
Il993al ). growing the NCO disk population (N+), particu- 
larly at small disk radii. Resonant relaxation and the Kozai 
mechanism will also allow non-disk NCOs to trade eccentric- 
ity with inclination and migrate into the disk over time (e.; 
iRauch fc Tremaindl 19961 ; ISubr fc Karasll2005l ; IChandl2009 
On the other hand, NCO interactions within the disk that 
lead to ejection should keep ejected NCOs at relatively low 
inclinations, thereby increasing the probability of re-capture 
by the disk (growing N+). A disk capture rate of ~ 10 -4 of 
the non-disk NCO population over the lifetime of the AGN 
disk, corresponds to the capture of ~ 1O 3 M0 of non-disk 
NCOs mostly in the inner AGN disk, over the ~ 50Myr disk 
lifetime around a 1O 8 M0 supermassive black hole. 

Unlike protoplanets in disks around stars, a large num- 
ber of NCOs in AGN disks should have retrograde orbits. 
The behaviour of retrograde NCOs will depend on the effi- 
ciency of angular momentum transfer between the NCO and 
the disk gas. On one hand, if the coupling between NCO and 
gas is strong, the retrograde NCO rapidly loses angular mo- 
mentum and falls into the central supermassive black hole 
very quickly. On the other hand, if the coupling is very weak, 
the disk gas can move fast enough past the NCO that the 
NCO can persist in the disk for a long time without migrat- 
ing. In this case, prograde NCOs will migrate and encounter 
retrograde NCOs. 



4.2 Type II IMBH migration 

An IMBH that grows large en ough by accreting g as can 
open a gap in the AGN disk (|Sver fc Clarke! 1 199a ) . This 
phenomenon is analagous to gap opening by massive plan- 
ets in protoplanetary disks (|Armitagell2010l ). For typical disk 
parameters (H/r ~ 0.05, a = 0.01), the mass ratio required 
to open a gap is q ~ 10 -4 . However, there is a strong depen- 
denc y on disk viscosity for both the profi le and depth of the 
gap i|Crida et alj|2006l ; iMuto et alj|2010h. To open a g ap in 
the disk requires low disk viscosity ( Crida et al.l l2006) 



a < 



(23) 



wher e a is the disk viscosity parameter ( Shakura & Sunyaev 
1973). In the AGN disk modelled by ISirko fe Goodmanl 
(2003), H/r ~ 0.05 on average between 10 2 — 10 5 r o and 



l|Sver fc Clarkelll995l ). Thus, a gap open ed by a 3 x 10 4 M B 
IMBH in the AGN disk modelled by ISirko fc Goodmanl 
(2003) will be closed by pressure and/or accretion in the 
outermost and innermost parts of the disk where (H/r) and 
a are large. In a more viscous type of accretion flow (a ^ 0.1, 
e.g. advection dominated), an IMBH might not open a gap 
in the disk. Whether an IMBH opens a gap will have ma- 
jor implications for observational signatures in AGN, but we 
defer that discussion to Paper II. 

An IMBH that opens a gap will tend to migrate on the 
viscous disk timescale (Type II migration) given by 



Til = 



-i-Y 1 - 

a \r I uj 



(24) 



a = 0.01 is fixed. An IMBH with q > 3x 10~ 4 (> 3x 10 4 M Q ) 
will clear a gap in this disk. The gap will close by pressure 
if (H/r) > (q/a) 1/2 and by accretion if (H/r) > (q 2 /a) 1/5 



From eqn. [24] for a ~ 0.01 and H/r ~ 0.05 (approxi- 
mately the condit i ons a cross the AGN disk modelled by 
ISirko fc Goodmanl (|2003l )). the Type II migration timescale 
is ~ 10 4 x the orbital timescale. So, at 10 4 (10 5 )r 9 , the 
Type II migration timescale is ~ l(30)Myrs. Evidently a 
gap-opening IMBH can migrate on timescales shorter than 
the AGN lifetime across the disk. We therefore have to ask 
whether any gap-opening IMBH will survive the AGN disk? 
This is analagous to a major problem encountered in proto- 
planetary disk theory. The migration of some gap-opening 
migrators must somehow stall before accretion onto the cen- 
tral mass, otherwise no Jupiter-mass planets would be ob- 
served. One solution to this problem is that Type II mi- 
gration can stall once the co-rotating disk mass is less than 
the migrator mass. This condition could arise due to disk 
drainage onto the supermassive black hol e, or a change in the 
surfac e density profile of the disk. In the lSirko fc Goodmanl 
(2003) disk, this radius is ~ 10 4 r 9 for a 3 x 1O 4 M IMBH. 
Once Type II migration stalls, it can resume but at a much 
slower rate, once the migrat ors' angular momen t um is ex- 
porte d to the local disk (e.g. ISver fc Cl arke 1995; lArmitagel 
l2010h . 

If we simply assume the IMBH can undergo Type II in- 
migration without stalling, the IMBH will collide with NCOs 
at radii interior to its starting position 10 4 (10 5 )r 9 before ac- 
creting onto the supermassive black hole in l(30Myrs). Ig- 
noring gas accretion, the IMBH can swallow up to 5%(50%) 
of the uniformly distributed disk NCO population if it starts 
migrating at 10 4 (10 5 r 9 ). So, IMBH growth via NCO merger 
can be as much as ~ 5OOOA/0 /30Myr (starting at 10 r g and 
10 4 NCOs in the disk). This growth rate is due to NCO 
mergers only and does not include growth due to gas ac- 
cretion (see H5. 21 below). If the IMBH stalls permanently at 
~ 10 4 r 9 , only a small number of remaining migrators (~ 2% 
of the remaining disk NCO population) migrate inwards to 
merge with the IMBH within 50Myrs. 

To sum up, for a powerlaw stellar mass function (~ 
M" 3 ), most NCOs should be low mass stars (< IMq) with 
a small fraction of compact NCOs (mostly white dwarfs, 
some neutron stars and stellar mass black holes) . The largest 
mass NCOs (and seeds for IMBHs) are likely to be small in 
number, and consist of stellar mass black holes or massive 
binaries. IMBH seedlings will undergo Type I migration in 
the disk. IMBH migration means that they can maintain a 
feeding zone of disk NCOs as they 'catch up' with the much 
more slowly migrating low mass disk NCOs. This migration 
of the feeding zone is precisely analagou s to the situation ex - 
pected for migrating giant planets (e.g. lAlibert et al"1l2004l ). 
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Some of the classic problems of protoplanetary migration 
(e.g. how to stop migrators from accreting onto the central 
object, or how to get migrators moving once stalled) will 
also apply to IMBH seeds in AGN disks. Nevertheless, low 
relative velocity encounters due to migration will result in 
large NCO collision cross-sections and will help grow IMBH 
seeds via core accretion. 



5 A MODEL OF IMBH GROWTH IN AGN 
DISKS 

In this section, we shall draw together much of the above 
discussion and construct a simple model of IMBH growth 
in AGN disks. The starting point for our model is a stel- 
lar mass black hole. This 1OM black hole can accrete gas 
from the AGN disk, migrate within the disk and collide 
with disk NCOs. We follow the approac h of models of giant 
planet growth in p rotoplanetary disks (|Pollack et al.|[l996l ; 
lAlibert et al J 120040 . In section [BTTl we discuss the growth of 
the IMBH seedling as a result of collision with disk NCOs 
within the IMBH feeding zone. In section IST21 we discuss the 
growth of the IMBH seedling as the result of gas accretion. 



5.1 The parallel with 'core accretion' 

We considered a simple model of the growth of a lOAf 
IMBH seed embedded in an AGN disk around a 1O 8 M 
black hole. The disk NCO initial population is (1O 3 )1O 4 M , 
with a mass function of dN/dM oc M~ s (as discussed above 
mostly 0.6-Mq stars). We assume that the IMBH seed 'feeds' 
on NCOs within its accretion zone, given by a ± Sa where 
8a ~ 4Rh, analagous to the feeding zone of giant planet 
cores (jPollack et al.l 0*9961 ) . The maximum (isolation) mass 
that the IMBH seed can attain by feeding on all the NCOs 
within its accretion zone is 



with Vex, ~<r~<e >' Vk, the rate of IMBH seed growth 
via core accretion within its feeding zone is 



Miso = Mi + 167to 2 E f |) 



1/3 



(25) 



where Eo is the mean initial NCO surface density, q is the 
mass ratio of the IMBH seed to the supermassive black hole. 
If there are 10 M© NCOs initially in the disk around a 
10 8 Mq black hole, the mean initial NCO surface density is 
E ~ 3.5g/cm 2 . For an IMBH seed of Mi = 1OM , eon. [231 
implies Miso ~ 10 + 9(900)M Q at 10 4 (10 5 )r 9 . So, in prin- 
ciple, a stellar mass black hole in the outer disk could grow 
to many times its original mass just by accreting low mass 
disk NCOs. This process is analagous to the growth of giant 
planet cores by planetesimal accretion (jPollack et alj|l996t 
lArmitagdl20lch . 

Assuming small eccentricities (< e 2 > 1/,2 = Aa/a) for 
disk NCOs, we are in a shear-dominated regime and the rate 
of mass growth of the IMBH seed may b e approximated by 
the form of giant planet core growth as (|Armitagell2010l ) 



dM _ 9 (Aaf 
~dt ~ 32 < i 2 >V2 a R H 



(26) 



where < i 2 > 1//2 is the rms inclination for the NCO distri- 
bution, v is the relative local overdensity of disk NCOs and 
Cecil is the collision cross-section as given by eqn.fj] Thus, if 
Aa =< e 2 > 1/2 a and if we assume < e 2 > 1/2 ~ 2 < i 2 > 1/2 , 



dM _ 9 vT, nTTr p 2GM 
IT ~ 8e(q/3y/3 v 2 



(27) 



where r v is the periastron. The periastron for compact 
object collisions with a n IMBH seed is r p ~ Rq 
(|Miller fc Hamilton! 120021 '). Substituting into eqn.~~7] where 
we assume < e 2 > 1/,2 ~ 0.01 is the equilibrium eccentricity, 
we find that for a 1OM IMBH seed at ~ 2 x 10 4 r 9 with 
E = 3.5gcm" 2 and v = 1, then dM/dt ~ lO~ 7 M /yr, 
which is approximately half the Eddington rate of growth. 
This is a very high accretion rate for a black hole, exceeding 
inferred accretion rates from gas disks in most Seyfert AGN 
l|McKernan et alj l2007h . Of course Eo is the average sur- 
face density assuming a uniform distribution of disk NCOs 
and that the mass ratio of disk NCOs to gas in the disk is 
~ 1%. If, instead the mass ratio is a factor of a few larger, 
or the surface density distribution of NCOs is non-uniform, 
the accretion rate of disk NCOs by IMBH can be substan- 
tially super-Eddvngton. Equally, if gas damping is more effi- 
cient than outlined above so that equilibrium is reached at 
e < 0.01 (e.g. due to a lower ratio of cusp population to 
disk NCOs), we could also reach super-Eddington rates of 
IMBH growth via mergers. From our earlier discussion, it 
is easy to envisage regions of the disk where NCOs tend to 
pile-up, leading to non-uniform distributions of disk NCOs. 
For example, if there is a region of the disk where inward 
an d outward torques balance (as in the scenario outlined 
bv lHorn et all (2012)), or migration stalls due to a change 
in the aspect ratio, or the mass of co-rotating gas drops. 
In these cases, we could write Eo ~ i/3.5gcm 2 , where v is 
an overdensity factor (whi ch could locally b e > 100 in a 
pile-up scenario such as in iHorn et al l l| 20121 ) and lead to 
highly super-Eddington growth). Note that a growing IMBH 
seedling avoids problems in protoplanetary coagulation the- 
ory, such as fracturing and sticking efficiency. For IMBH 
seedlings , nearby objects will either be captured (at high 
efficiency) or they will escape. 

Figure [4] shows the analytic growth of a 1OM stellar 
mass black hole (IMBH seed) and a 100M Q IMBH in an 
AGN disk around a 10 s M supermassive black hole. The 
lower curve in each case corresponds to v = 1 (no overden- 
sity, fiducial numbers) and the upper curves corresponds to a 
moderate over-density v ~ 5 of disk NCOs (or equivalently, 
a slight overdensity, v — 2, and a moderately lower mean 
eccentricity e = 0.004). The fiducial mass doubling time for 
a black hole accreting gas at the Eddington rate (assuming 
10% efficiency) is 4 x 10 7 yrs. In the case of IMBH seeds grow- 
ing via collisions in AGN disks, and assuming a gas accretion 
rate of ~ Eddington, the total growth rates are x 1.5(3.5) 
Eddington for v — 1(5). At 3.5x Eddington growth rates, 
the mass doubling time could be as little as ~ llMyrs. Once 
the IMBH reaches its isolation mass, it will then grow via 
accretion from the gas disk, at a much slower rate (Edding- 
ton or a fraction thereof). However, with an increased mass, 
the IMBH will have a shorter Type I migration timescale 
(see ig] above). Thus, in ~ llMyrs, the IMBH will have mi- 
grated inwards or outwards in the disk and the size of the 
feeding zone will have grown. So, a super-Eddington mode 
of accretion via collisions could continue (dashed curve in 
Fig. ~4j . This process is analagous to the migration of giant 



10 B. McKernan, K.E.S. Ford, W.Lyra & H.B. Perets 



planet cores and th eir feeding zones in protoplanetary disks 
l|Alibert et alj|2004h . 

5.2 IMBH and gas accretion: runaway growth 

A big difference between IMBH growth in stellar clusters 
and in AGN disks is that in the latter, gas can damp orbits 
quite effectively. So mergers tend to be more frequent, and 
the IMBH seeds can continuously accrete dense gas. Thus, 
we expect IMBH growth in AGN disks to involve growth 
by merger (as in stellar clusters) but we also expect growth 
by gas accretion. Torques from the gas will cause the IMBH 
to migrate and enter new feedin g zones, analagous to the 
situation in protoplanetary disks (|Alibert et al.ll2004r i. If the 
IMBH grows large enough (q = 1(T 4 or 1O 4 M around a 
1O 8 M0 supermassive black hole), the IMBH can open a gap 
in the gas disk and the rate of gas accretion will drop. Note 
that although we concentrate on building a single IMBH, 
multiple IMBH seedlings (10 — lOOM©) are likely to appear 
in the disk (assuming dM/dt oc M -3 , see discussion above). 

One problem will be in preventing an IMBH from mi- 
grating onto the supermassive black hole. Outward migra- 
tion and the stalling of migration due to a drop in disk 
surface density or a change in the disk aspect ratio are pos- 
sible solutions, but as with protoplanetary disk theory, this 
theoretical problem is complicated and remains unsolved for 
now. A sufficiently massive gap-opening IMBH O 10 4 A/ Q 
around a 1O 8 M0 supermassive black hole) will grow if the 
AGN disk is particularly long-lived (> 50Myrs), or if there 
is a large local disk NCO overdensity (y), or if gas damp- 
ing is particularly efficient so that equilibrium eccentric- 
ity is e < 0.01. Alternatively, an IMBH which survives 
a period of AGN activity could grow to gap-opening size 
via the mechanisms outlined here, during a later, indepen- 
dent AGN phase. Earlier in the history of the Universe (at 
z ~ 2), the time b etween individual AGN phases should be 
much smaller (e.g. lyittorini et ai1l2005l; IPohertv et al-lboOrj ; 
IShankar et"ail2 009: McKernan ct al. 2 oTobT) . So, if they sur- 
vive, large mass (gap-opening) IMBH could grow rapidly in 
galactic nuclei over a few 100 Myrs and observational signa- 
tures of IMBH in galactic nuclei may be common at higher 
redshift (see Paper II). 

Of course, as the gas disk is consumed or blown away, 
other mechanisms will come into play. For planetesimals in 
a late-stage protoplanetary disk, planet-planet interactions, 
the Kozai mechanism or resonant relaxation can increase e, i 
(e.g. iRauch fc Tremaine! Il996l : ISubr fc Karasl 120051 : IChand 
2009). By analogy with protoplanetary disks, such mecha- 
nisms will certainly apply in the late stages of an AGN disk 
when most gas has been drained. However, we do not con- 
sider these mechanisms in more detail here since damping 
due to dense gas disk should dominate such effects (see Pa- 
per III for further discussion). 



6 CONCLUSIONS 

We show that it is possible to efficiently grow intermediate 
mass black holes (IMBH) from stars and compact objects 
within an AGN disk. Nuclear cluster objects (NCOs) in the 
AGN disk are subject to two competing effects: orbital exci- 
tation due to cusp dynamical heating and orbital damping 
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Figure 4. The growth over time of a IOMq IMBH seed and 
a 100M Q IMBH located 2 X 10 4 r s from a 10 8 M Q supermassive 
black hole in an AGN disk. The solid (dashed) curves show IMBH 
growth as it accretes disk NCOs within ^ Ar^ at a rate x0.5(2.5) 
Eddington, assuming e = 0.01 and an initial NCO surface den- 
sity So ~ 3.5(18)gcm -2 . The dashed curves could equivalently 
be generated with Sq ~ 7gcm~ 2 and e = 0.004. The gas ac- 
cretion rate is assumed to be at the Eddington rate, for a total 
accretion rate of X 1.5(3.5) Eddington and a mass doubling time 
of ~ 27(ll)Myr respectively. We assume the IMBH migrates in 
the disk and moves its feeding zone, so that it may continue to 
undergo collisions with NCOs at up to supcr-Eddington rates 
(anal agous to continued c ollision and migration in protoplanetary 
disks lAlibert et al.l 1120041) 1. In order to clear out > 1O 3 M0 of disk 
NCOs from this choice of initial disk location, outward migration 
is required. Once the IMBH reaches 1O 4 M0 we assume the disk 
NCO population has been cleared out and the IMBH opens a gap 
in the disk, accreting gas at an Eddington rate thereafter. In this 
picture, between 1/3 and 2/3 of the IMBH mass is actually due 
to gas accretion. 



due to gas in the disk. For a simple, semi-analytic model 
we show that gas damping dominates such that equilibrium 
eccentricities of disk NCOs are e ~ 0.01. In this case IMBH 
seedling formation via NCO collision is more efficient in the 
AGN disk than in stellar clusters (the standard model for 
IMBH formation). If, as we expect, gas damping dominates 
then the dynamical heating of disk NCOs by cusp stars is 
transmitted to the gas disk. This is a new, additional source 
of heating of the outer disk that can help counter the well- 
known gravitational instability (Q ^ 1) of the outer disk. 

Stellar mass black holes and hard massive binaries are 
likely IMBH seeds. IMBH seedlings grow by collisions with 
disk NCOs within their feeding zone (a ± 4Rh) at near Ed- 
dington rates, as well as via gas accretion. IMBH seedlings 
will migrate within the AGN disk and so continue to feed on 
disk NCOs as they accrete gas. If there are regions of mod- 
est over-density of NCOs in the disk, IMBH seedling growth 
via collisions can be super-Eddington and a 10M Q IMBH 
seed orbiting a 1O 8 M0 supermassive black hole can grow to 
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~ 3OOM0 in less than the fiducial AGN disk lifetime. An 
over-density of disk NCOs can occur in regions of the disk 
where e.g. outward torques and inward torques balance, or 
where the aspect ratio changes, or where IMBH migration 
stalls. 

The largest IMBH will open gaps in AGN disks, 
analagous to giant planets in protoplanetary disks. Gap- 
opening IMBH are more likely to arise if: gas damping is very 
efficient (equilibrium disk NCO eccentricity is e < 0.01), or if 
the disk is long-lived (> 50Myrs), or disk NCO surface den- 
sity is moderately high (> 15gcm~ 2 ), or if there is an IMBH 
seedling which survived a previous AGN phase (analagous to 
the survival of planets in protoplanetary disks). Our model 
of IMBH growth in AGN disks strongly parallels the growth 
of giant planets in protoplanetary disks. We leave a discus- 
sion of model predictions, observational constraints and im- 
plications of efficient IMBH growth in AGN disks to Paper 
II. 
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